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We investigate the two-dimensional magnetic Schrodinger operator 

^ ■ Hb,i3 = (— iV — A)"^ — P6{- — r), where F is a smooth loop and the 

. vector potential A corresponds to a homogeneous magnetic field B 

^ , perpendicular to the plane. The asymptotics of negative eigenvalues 

^ I of H^^fs for /? — > oo is found. It shows, in particular, that for large 

I enough positive /3 the system exhibits persistent currents. 

si 

1 Introduction 

One of the most often studied features of mesoscopic systems are the per- 
sistent currents in rings threaded by a magnetic flux - let us mention, e.g., 
||(JG1^, |(JWB|] and scores of other theoretical and experimental papers where 



they were discussed. For a charged particle (an electron) confined to a loop F 
the effect is manifested by the dependence of the corresponding eigenvalues 
Xn on the flux threading the loop, conventionally measured in the units of 
flux quanta, 27rhc\e\~^. The derivative dXn/dcf) equals —-In, where is the 
persistent current in the ra-th state. In particular, if the particle motion on 
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the loop is free, we have 



(21,-'' 



where L is the loop circumference and m* is the effective mass of the electron, 
so the currents depend linearly on the applied field in this case. 

Of course, the above example is idealized assuming that the particle is 
strictly confined to the loop. In reality boundaries of a quantum wire are po- 
tential jumps at interfaces of different materials. As a consequence, electrons 
can be found outside the loop, even if not too far when we consider energies 
at which the exterior is a classically forbidden region. 

A reasonable model respecting the essentially one-dimensional nature of 
quantum wires is a 2D Schrodinger operator with an attractive 5-interaction 
on an appropriate curve F, or more generally, a planar graph. Since the 
interaction support has codimension one, the Hamiltonian can be defined 
through its quadratic form and the corresponding resolvent can be written 
explicitly as a generalization of the Birman-Schwinger theory ||B'1|, |BEKIS 



This leads to some interesting consequences such as the existence of bound 



states due to bending of an infinite and asymptotically straight curve 



A natural question is how such a model is related to the ideal one in which 



the electron is strictly confined to the curve F. In [ [EY|| we have derived an 
asymptotic formula showing that if the b coupling is strong, the negative 
eigenvalues approach those of the ideal model in which the geometry of F is 
taken into account by means of an effective curvature-induced potential. The 
purpose of this paper is to ask a similar question in the situation when the 
electron is subject in addition to a homogeneous magnetic field B perpendic- 
ular to the plane. We are going to derive an analogous asymptotic formula 
where now the presence of the magnetic field is taken into account via the 
boundary conditions specifying the domain of the comparison operator. 

An easy consequence of this result is that for a strong enough (5-interaction 
the negative eigenvalues of our Hamiltonian are not constant as functions of 
5, i.e. that the system exhibits persistent currents. Their further properties 
depend, of course, on the specific shape of F; this fact as well as the stability 
of such currents with respect to a disorder raise questions about optimal ways 
of interpreting the corresponding magnetic transport. We comment on this 
point in the concluding remarks. 
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2 Description of the model and the results 



As we have said above we are going to study the Schrodinger operator in 
L^(]R^) with a constant magnetic field and an attractive 5-interaction on a 
loop. For the sake of simplicity we emply rational units, h — c — 2m* — 1 
and absorb the electron charge into the field intensity B. We shall use the 
circular gauge, A{x^y) = [—^By,^Bxy 

Let r : [0, L] 3 s ^ (ri(s), r2(s)) G be a closed counter-clockwise 
Jordan curve which is parametrized by its arc length. Given /3 > and 
B eR,we define 



(^-id, + ^By^f + (^-idy-^Bx^f -pj^ 



'ds 



with the domain i/^(R^), where = d/dx etc., and the norm refers to 
L^(]R^). It is straightforward to check that the form qs/s is closed and below 
bounded. We denote by Hb^p the self-adjoint operator associated with it 
which can be formally written as 



H 



B,I3 



-iV -Af -i35{--T). 



Our main aim is to study the asymptotic behavior of the negative eigenvalues 

of i?B,/3 as /5 — s> +00. 

Let 7 : [0, L] 9 s ^^ (r'/r'2 - r'2T;)(s) G M be the signed curvature of T. 
We denote by VL the region enclosed by F, with the area \Vt\, and define the 
operator 

Sb--^, - -Ms? 



on L^((0, L)) with the domain 



where y?*-'^'-' stands for the A'-th derivative. 

We fix j e N and denote by fJ^j{B) the j-th eigenvalue of Sb counted with 
multiplicity. Our main results is the following claim. 

Theorem 2.1 Let n be an arbitrary integer and let ^ I G M. be a compact 
interval Then there exists P{n,I) > such that 



U<^d{HB,(3)r\ {-00,0)} >n for (5 > (5{nJ) and Bel. 
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For (3 > I3{n) and B E I we denote by Xn{B,(3) the n-th eigenvalue of Hs^p 
counted with multiplicity. Then Xn{B,f3) admits an asymptotic expansion of 
the form 

XniB,l3) = -^(3'' + finiB) + Oil3-Hnl3) as /3 ^ +00 , 
where the error term is uniform with respect to B E I. 

Recall that the flux (p through the loop is i?|i7|/27r in our units. The existence 
of persistent currents is then given by the following consequence of the above 
result. 

Corollary 2.2 Let ^ ^ I G be a compact interval and let n E N. Then 
there exists a constant i3i{n,I) > such that the function A„(-,/5) is not 
constant for [3 > Pi{n, I). 



3 The proofs 



Since the spectral properties of Hb,i3 are clearly invariant with respect to 
Euclidean transformation of the plane, we may assume without any loss of 
generality that the curve F parametrizes in the following way, 

PS PS 

ri(s) = ri(0) + / cosH{t)dt, r2(s) = 12(0) + / smH{t)dt, 
Jo Jo 

where H{t) := — "y{u) du. Let \E'a be the map 

: [0,L) X (-a,a) 3 {s,u) ^ iTi{s)-uT'2{s),T2{s)+uT[{s)) E Ml 



By |[EY| , Lemma 2.1] we know that there exists ai > such that the map 
is injective for all a E (0,ai]. We fix thus a E (0,ai) and denote by the 
strip of width 2a enclosing F 

:= ^a{[0,L) X (-a, a)) . 

Then the set ]R^\Sa consists of two connected components which we denote 
by A™ and A"*^*, where the interior one. A™, is compact. We define a pair of 
quadratic forms, 

2 



aAf] 
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(^-td, + ^By^ f 


+ 







-idy - -Bx f 



-P j^\f{x)\'ds.. 



which are given by the same expression but differ by their domains; the latter 
is i?o(E„) for and i^H^a) for 

quadratic forms 



Furthermore, we introduce the 
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for j = out, in, with the domains if|(A{) and H^{Al) corresponding to the ± 
sign, respectively. Let ^ ^'b^o. i ^^"^ -^B,'a t>e the self-adjoint operators 
associated with the forms ^ ^b^I'^' ^^"^ ^s'a ' respectively. 

As in | [li)Y|| we are going to use the Dirichlet- Neumann bracketing with ad- 
ditional boundary conditions at the boundary of S^. It works in the magnetic 
case too as one can see easily comparing the form domains of the involved 



operators - cf. Thm. XIII.2]. We get 



E 



B,a 



E 



B,a 



E 



out,- 
B,a 



(3.1) 



^B,a ^ B,a,f3 

with the decomposed estimating operators in L^(]R^) = L^(A™) © L^(Sa) © 
L^(A°"*). In order to assess the negative eigenvalues of Hb,p, it suffices to 
consider those of L^^^p ^^'^ ^~Bap^ because the other operators involved in 
( |3.1| ) are positive. Since the loop is smooth, we can pass inside to the 
natural curvilinear coordinates: we put 

{Uaf){s,u) = {l + U^{s)Y/^f{^a{s,u)) for f E {Jla) 

which defines the unitary operator Ua from L^(Ila) to L^((0,L) x (—a, a)). 
To express the estimating operators in the new variables, we introduce 

Ql = {ipeH\{Q,L)x{-a,a))- ip{L, ■) = ip{Q, ■) on (-a, a), 

¥?(-,a) = (/?(-,-a) = on (0,L)|, 
Q- = ^ipeH\{Q,L)x{-a,a))] ip{L, ■) = ip{Q, ■) on (-a, a) 

and define the quadratic forms 



B,a,l3 



[g] 



(3.2) 



(1+-U7(s)) \dsg\ duds + 



\dug\^ duds 
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-L 

V{s,u)\g\^dsdu- 13 I \g{s,0)\'^ds 
Jo 

b± /"^ lis) , , ,,2^ 6± /-^ 7(s) , , ,,2^ 

+ ^ /" /" B\Tj-2uT^r^ + Tl + 2uT2T[ + u^)\g\'^duds 

+- 



a 
L 



(T2 + uT[) ((1 + M7)-^ cosHgdsg - sin H gdug) duds 

-a 

(Fi — Mr2) ((1 + M7)~^ sinHgdsg + cos H 'gdug) duds 

a 

on 5 respectively, where 6+ = and 6_ = 1, and 

V^(.,«) = ^(l+«7(^))-W'(^)-^(l+«7(s))-Vy(s)^-i(l+ti7(^))-'7(^)' 



is the well-known curvature- induced effective potential \^S\\ . Let D"^ a /3 t>e 
the self-adjoint operators associated with the forms b^^^, respectively. In 
analogy with |[EY| , Lemma 2.2], we get the following result. 



Lemma 3.1 U*D%^^ JJa = 1%^^^. 

The presence of the magnetic field gave rise to terms containing gdgg and 
'gdug in ( p.2|) . In order to eliminate the corresponding coefficients modulo 
small errors, we employ another unitary operator. We put 



^-B I\r2{t)r[{t) -r'^imit)) dt- 

it follows from the Green theorem that Tb{L) = B\Q\ . Then we define 



Tsis)-. 2- 



r 1 ~1 

(MBh){s,u) := exp iTs (s) + -Bu {r2{s) sin H{s) + ri{s) cos H{s)) h{s,u) 

for any h G L'^{{0, L) x (—a, a)); it is straightforward to check that Mb is a 
unitary operator on L^((0,L) x (—a, a)). We define 

Qta = {v^H\{Q,L)x (-a, a)); <^(L, ■) = e-*^l^l<^(0, ■) on (-a, a), 

V9(-,a) = V9(-, -a) = on (0,L)|, 
QB,a = {^Gi7i((0,L)x(-a,a)); ^.(L, ■) = e-^l^l¥;(0, •) on (-a, a)}. 



and another pair of quadratic forms 

+ [B{r2 + ur[){l + u-f)-^ cosH - B(ri - 'ur^)(l + -^7)-^ sin H 
-B{1 + u-f)-^(r2 COS H-Ti sin H) 

+B{1 + u^)-^{r2smH + Ti cos HYu] lm{gdsg) + WB{s,u)\g\^^ duds 

for g G a, respectively, where 
Wb{s,u) 

= V{s, u) + ^il + ii7)-2S V((r2 sin H + r^ cos H)')^ 

+^B\Tl - 2uTiT'2 + Tl + 2Mr2r; + u^) 

+B(r2+ur[)(l + u-f)-^T^(s)cosH - B(ri-ur'2)(l + u-f)-^T'B(s) sin// 
+^(1 + u-f)-'^B'^{r2 cos H-Fi sin H)'^ + ^S^(r2 sin H + Fi cos //)^ 

+ [B{T2 + uT[){l + u-f)-\osH - B{T I- uT 2) (l+u-f)-^ sin H 
-B{1 + U'y)-\r2 cos H-Fi sin //)] -5(r2 sin // + Ti cos HYu 

+ [-5(r2 + -ur;) sin // - S(ri - 'wr'2) cos //] ^S(r2 sin // + Ti COS //) . 

Let Dsap tie the self-adjoint operators associated with the forms ^fa^, 
respectively. By a straightforward computation, one can check the following 
claim. 

Lemma 3.2 M* /^± „_^Mb = D± . 

The next step is to estimate /^^a/3 operators with separated variables. 
Denoting 

7+: = maf|7(-)| 



we put 



Nsia) := max 

{s,u)e[0,L]x[-a,a] 



B{T2 + uT[){l + u'y)-^cosH 

-S(ri - -ur'J (1 + -^7)-^ sin 
-B{1 + u-f)-^{r2 cos H-Fi sin H) 

+B{1 + U7)"^(r2 sin // + Ti cos Hyu 



and 



Msia) :— max 

{s,u)e[0,L]x[-a,a] 



WB{s,u) + -^{sf 



Let 7^ / C M be a compact interval. Then there is a positive K such that 



Neia) + Meia) < Ka for < a < 



1 



27+ 



and Bel, 



where K is independent of a and 5. For a fixed < a < 3:^^, we define 




J -a 



[l^a^^y±-NB{a) 



+ 



/? r\f{s,0)\'ds-^+b^ r{\f{s,a)\'+\f{s,-a)\')ds 
Jo Jo 



-^j{sf±^NB{a)±MBia) 



for / e Qsa' respectively. Since \lm{gdsg)\ < l^sS'l^), we obtain 



^^,a,,[/] < bt^ajf] for /GQ^,,, 



(3.3) 

(3.4) 



Let Hsap ^® self-adjoint operators associated with the forms &^a/3' 
respectively. Furthermore, let T^^ be the self-adjoint operator associated 
with the form 



\f'{u)\'ds-(3\f{0)\\ feH',{{-a,a)), 
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and similarly, let g be the self-adjoint operator associated with the form 



We define 



\f\u)\'ds-P\m\'-M\m\'+\f(-a)n / e H\(-a,a)) 



U 



B,a 



:iTa-f+y'±-NBia) 



^-l^(sy±^NB{a)±MB{a) 



in L^((0, L)) with the domain Pb specified in the previous section. Then we 
have 

^l,a,/J = t^l,a®l + l®7t/9- (3-5) 

Let 11^(3, a) be the j-th eigenvalue of C/^^ counted with multiplicity. We 
shall prove the following estimate. 

Proposition 3.3 Let j e N. Then there exists C{j) > such that 
\l^+{B,a)-i^j{B)\ + \i^j{B,a)-i^j{B)\ < C{j)a 



holds for B & I and < a < where C{j) is independent of B and a. 
Proof: Since 

1 



1 - 07+)-' + -Nsia) 



a7+(2 - a7+) 1 
(l-a7+)2 +2 



NB(a) 



and since A^B(a) + Msia) < Ka for < a < and B e I, we infer that 
there is a constant Ci > such that 



[1 - aj+)-^ + -Nsia) 



for < a < 2^ and B & I. This together with the min-max principle implies 
that 

/x;(5,a)- (l-a7+)-2 + liVB(a) iXj{B) < Cia 
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for < a < 2^ and B E I. Since is continuous, we claim that there 
exists a constant C2 > such that 

\fij{B,a)-fij{B)\ < C2a 

for < a < 2^ and B E I. In a similar way, we infer the existence of a 
constant C3 > such that 

\fij{B,a) - fij{B)\ < Csa 

for < a < and B e I. M 

27+ 

We also recall the following result from [EY]. 

Proposition 3.4 (a) Suppose that /5a > |. Then T^^ has only one negative 
eigenvalue, which we denote by (^p. It satisfies the inequalities 

< < + 2/5' exp (-i/3a) . 

(b) Let aP > 8 and (3 > |7+. Then has a unique negative eigenvalue 
C,~p, and moreover, we have 



± 



Now we are ready to prove Theorem prT| . We put a{[3) = 6/? In /5. Let ■ 

(/3),/3 



be the j-th eigenvalue of ^^^^ ^, by Proposition we have 

"it /■± f± 



/^From the decompositions ( p.5|) we infer that {^^ j + «(/5))}j,A;GN5 prop- 

erly ordered, is the sequence of the eigenvalues of ^^^-^ ^ counted with 
multiplicity. Proposition p]^ gives 



ej, + am > /if (S, a(/5)) = ^l^iB) + 0{p-' In /3) (3.6) 

for B E I , j > 2, and /c > 1, where the error term is uniform with respect to 
B E I. For a fixed j e N, we put 
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Combining Propositions O and O we get 

rlp^, = -\(3^ + H{B) + 0{(3-Hn(3) as /3 ^ oo , (3.7) 

where the error term is uniform with respect to -B G /. Let us fix now n G N. 
Combining ( |3.6| ) with ( p.7| ) we infer that there exists (3{n,I) > such that 
the inequahties 

< 0' ^B,l3,n < ^tj + f^ki^, a(/5)), r^^^ „ < + /i^ (5, a(/?)) 

hold for B & I, j3 > l3{n,I), j > 2, and > 1. Hence the j-th eigenvalue 
of -^^a(^)^ counted with multiplicity is 13 j for B ^ I, j < n, and /? > 
/?(n,/). Let 5 G / and /3 > (3{nJ). We denote by K,f{B,(3) the j-th 
eigenvalue of L^^b- Combining our basic estimate with Lemmas |3]I 



and relations ( |3.3| ) and (^.41), and the min-max principle, we arrive at 
the inequalities 

TB,(,,,<^^JiB,P) and <r+^^^. for l<j<n, (3.8) 

so we have k^{B,P) < < inf aess{HB,f3)- Hence the min-max principle and 
(|3.1| ) imply that Hb,i3 has at least n eigenvalues in (— 00, k+(_B, /5)]. Given 
1 < j < we denote by \j{B,P) the j-th eigenvalue of Hb^p- It satisfies 



Kj{B,[3)<X,{B,f3)<K+{B,[3) for l<j< 
this together with (|3.7|) and ( |3.8|) implies that 



n : 



Aj(5,/5) = -i/32 + /Xj(5) + C(/3~Mn/3) as /3 ^ 00 for l<j<n, 

where the error term is uniform with respect to B E I. This completes the 
proof. I 

Proof of Corollary ]^. 2^ By [[RS| , Thm XIII.89] no eigenvalue //„(■) is constant 
on /. This together with Theorem yields the claim. | 



4 Concluding remarks 

The above corollary answers the question we posed in the introduction as a 
mathematical problem showing that a ring with a strong enough attractive 
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^-interaction does exhibit persistent currents. On the other hand, from the 
physical point of view it would be bold to identify a mere non-constantness 
of the eigenvalues with a genuine magnetic transport around the loop. 

The problem is similar to other situation where an electron can be trans- 
ported in a magnetic field due to the presence of a "guiding" perturbation. 
A prime example are the edge currents ||Ha|, |MS|| which attracted a wave 



of mathematical interest recently in connection with the problem of stabil- 
ity of the transport with respect to perturbations. In case of a single edge 
and a weak disorder a part of the absolutely continuous spectrum survives 
||BP| , [b'GWI , [MMP|| but the fact itself gives no quantitative information about 



the transport. On the other hand, a system with more than one edge may 
have no continuous spectrum at all and still it has states in which electrons 



travel distances much larger than the corresponding cyclotron radius ||FM . 

In our case it is clear, for instance, that the loop geometry influences the 
transport substantially. If F is a circle, e.g., than up the C(/9~^ln/9) error 
the persistent-current plot will have the ideal saw-tooth shape as we can see 
from the relation ( |1.1| ); one expects that the eigenfunctions will be "spread" 
around the whole circle. In contrast, if the loop is rather "wiggly" the one- 
dimensional comparison operator Sb contains an irregular effective potential 
coming from the rapidly varying curvature, which may cause - depending 
on the strength of such a "disorder" - that the most part of the electron 
wavefunction will be concentrated in (the vicinity of) a part of the loop only. 
The same may happen if the loop curvature is slowly changing but a disorder 
potential is added to the Hamiltonian. 

To distinguish the situations with a significant tiSLnsport, one needs clear- 
ly to understand better the sketched "disordered" cases which do not fall into 
this category. We leave this problem to a future publication. 
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